Current reversals in a rocking ratchet: the frequency domain 
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Motivated by recent work [D. Cubero et al., Phys. Rev. E 82, 041116 (2010)], we examine 
the mechanisms which determine current reversals in rocking ratchets as observed by varying the 
frequency of the drive. We found that a class of these current reversals in the frequency domain are 
precisely determined by dissipation-induced symmetry breaking. Our experimental and theoretical 
work thus extends and generalizes the previously identified relationship between dynamical and 
symmetry-breaking mechanisms in the generation of current reversals. 
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I. INTRODUCTION 

In out-of-equilibrium systems, directed transport can 
be obtained without the application of a net bias force. 
Such a counter-intuitive phenomenon is usually termed 
the ratchet effect [iHl]. The occurrence of a directed 
current in such systems can be precisely related to the 
breaking of the relevant spatio-temporal symmetries 

El. 

An intriguing feature of ratchets are current reversals, 
where the sign of the generated current changes following 
a variation in a system parameter a spe- 

cific example, consider the case of a rocking ratchet, con- 
sisting of Brownian particles in an asymmetric sawtooth 
potential. A sinusoidal rocking force drives the system 
out of equilibrium, thus allowing the generation of di- 
rected transport. In such a system, current reversals can 
be observed by varying the noise strength at fixed ampli- 
tude of the rocking force, as well as by varying the force 
amplitude while keeping constant the noise strength [3]. 
The key feature of current reversals is that the sign of the 
current can be reversed by varying a system parameter, 
although the considered parameter does not change the 
symmetry of the Hamiltonian. At first sight this suggests 
that current reversals are a dynamical phenomenon, not 
traceble back to a symmetry-breaking mechanism. 

Previous theoretical work aimed to identify the mech- 
anisms underlying current reversals in rocking ratchets. 
For an underdamped deterministic ratchet, Mateos 13 1 
argued that current reversals induced by a variation of 
the driving strength correspond to a bifurcation from 
a chaotic to a periodic regime. However, the proposed 
mechanisms turned out not to be general, as in the very 
same system current reversals can be observed also in the 
absence of such bifurcations, and moreover not all chaos- 
to-order transitions necessarily lead to current reversals 
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A different approach toward the understanding of cur- 
rent reversals induced by a variation of the rocking force 
amplitude was recently introduced for a biharmonically 



driven spatially symmetric rocking ratchet [la l. It was 
shown that a class of current reversals is precisely deter- 
mined by symmetry breaking. In this way, a link was 
established between dynamical and symmetry-breaking 
mechanisms. The still open issue now is to which extent 
the above link can be generalized, either to different sys- 
tems or to current reversals of different type in the same 
set-up. 



In the present work, we generalise the link between dy- 
namical and symmetry-breaking mechanisms introduced 
in the aforementioned work. We consider the same 
ratchet system and examine the current reversals induced 
by a variation of the frequency of the rocking force. We 
notice that a variation of the rocking force frequency can- 
not be mapped onto a variation of the force amplitude, 
i.e. the current reversals considered here are not a pri- 
ori equivalent to the ones observed by varying the force 
strength. Our experimental and theoretical work shows 
that also in this case there is a class of current reversals 
that are determined by symmetry breaking, thus gener- 
alizing the argument put forward in Ref. ^161] . 



This work is organized as follows. In Section II we 
introduce the ratchet set-up used in this work, and re- 
call the elements of the symmetry analysis that will be 
needed to establish the link between current reversals in 
the frequency domain and dissipation-induced symmetry 
breaking. Section III presents experimental results ob- 
tained with cold atoms in a driven optical lattice in the 
regime of weak damping. Section IV analyses theoreti- 
cally the relationship between dynamical and symmetry- 
breaking mechanisms for the current reversals induced 
by a variation in the driving frequency. Conclusions are 
drawn in Sec. IVl 
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II. SET-UP AND SYMMETRIES 
A. Ratchet set-up 

In this work we consider a ratchet set-up consisting of 
Brownian particles in a spatially symmetric periodic po- 
tential driven by a time-asymmetric force. To capture the 
main features of the dynamics, in the theoretical analysis 
we consider the simple case of a linear friction, as usually 
considered in the ratchet literature. This will also allow 
us to explore different regimes of damping and noise. The 
relevant Langevin equation is in this case: 

mi = -ax-U'{x)+F{t)+^{t) (1) 

where U{x) — C/q cos(2fca;)/2 is a periodic potential, a 
is the friction coefficient, £,{t) is a Gaussian white noise: 
{£_{t)) = 0, (e(i)C(i')) = 21>(5(t - t'), and F{t) is a bihar- 
monic drive described by: 

F{t) = Fq [Ai cos{Lut) + A2 cos(2tjt + 4))] . (2) 

The appearance of a ratchet-effect, i.e. the generation 
of directed motion, in such a set-up is a very well estab- 
lished fact [iMl- 

B. Symmetry analysis 

We now recall the essential elements of the symmetry 
analysis that, initially developed to explain the genera- 
tion of a current [9-11], will allow us to introduce the 
basic concepts that will be needed to establish the link 
between current reversals in the frequency domain and 
dissipation-induced symmetry breaking. 

In general, the symmetry analysis [gj-fTTj is used to 
identify the symmetry of the system which prevent di- 
rected motion. For the specific case of a spatially sym- 
metric potential, of interest here, there are two of those 
symmetries: the shift symmetry, which corresponds to in- 
variance under the transformation (x,p, t) {~x, ~p, t+ 
T/2), with T the period of the drive, and the time- 
reversal symmetry, which requires invariance under the 
transformation {x,p,t) — > (x, — p, — t). A bi-harmonic 
drive of the form of Eq. [5] breaks the shift symmetry 
for any value of the relative phase 0. For the time- 
reversal symmetry it is necessary to distinguish differ- 
ent cases, corresponding to different levels of dissipa- 
tion. For no dissipation (Hamiltonian case), the sys- 
tem is symmetric under time-reversal for (/) — nn with 
n integer, and therefore for these values directed motion 
cannot be produced. It can be shown [2^ that the 
dependence of the average velocity on the phase (j) is, 
in leading order, v = Asiii(f). Consider now the case 
of non zero dissipation. For non zero dissipation, the 
time-reversal symmetry is broken by dissipation also for 
(f) = nn with n integer. Thus, a current can be gen- 
erated also for these values of the phase (j). For weak 
dissipation, the dependence of the average velocity on 



the phase (j) is, in leading order, 1; = Asm((j) — (po), where 
00 is a dissipation-induced symmetry-breaking phase lag 
which vanishes in the Hamiltonian limit [l^, [H, [l^l . Fi- 
nally, in the overdamped regime, the system is invari- 
ant under the so-called " supersymmetry" [ll| {x,p,t) ^■ 
{x + \/2, ~p, ~t), with A the spatial period of the poten- 
tial, for (j) = n/2 + nn, with n integer. For these values 
of the phase (/>, direction motion is not allowed. 



III. EXPERIMENTAL RESULTS 

The experimental set-up and procedure are substan- 
tially the same as the ones used in our previous work 
of Ref. and we recall here only the essential ele- 

ments. Up to 10® ®^Rb atoms are trapped and cooled 
down to ^ 50 fiK in a magneto-optical trap. The atoms 
are then loaded in a ID lin_Llin dissipative optical lat- 
tice, in which the atom-light interaction determines both 
a periodic potential for the atoms and the dissipation 
mechanism which leads to a friction force and to fluctu- 
ations in the atomic dynamics. A driving of the form of 
Eq. ^ is applied by phase modulating one of the lattice 
beams j25j. For all the measurements presented in this 
work, the ratio between harmonics is kept fixed: Ai = 2, 
A2 = 1. The motion of the atoms in the driven lattice is 
studied by imaging the atomic cloud with a CCD camera. 
The velocity of the center of mass of the atomic cloud is 
then derived by using these images. 

We first consider the standard configuration for the de- 
tection of current reversals: for a fixed Hamiltonian, with 
broken symmetry, the current is studied as a function of 
a parameter, whose variation does not change the sym- 
metry of the Hamiltonian. In the present case, we fix the 
ralative phase between driving harmonics </) = 7r/2, so to 
break the time reversal symmetry, and study the current 
as a function of the driving frequency uj. Results of our 
experiment are shown in Fig. [TJ The familiar situation 
of current reversals in the frequency domain 21,, 26] is 
observed: by varying the frequency of the drive it is pos- 
sible to revert the direction of the atomic current through 
the lattice. 

We now study the atomic current as a function of the 
phase (j) for different values of the driving frequency, in a 
range of frequencies around the value at which the cur- 
rent reversal is observed. This will allow us to establish 
a relationship between the observed current reversal and 
dissipative effects. 

The results of our measurements are shown in Fig. [21 
for a driving force amplitude corresponding to that of 
Fig- Wic). For all considered driving frequencies, the 
dependence of the current on the phase (j) is well described 
by v/vr — Asiii{(j3 — (/jq), with 0o a dissipation-induced 
symmetry-breaking phase l ag. This is in agreement with 
previous observations [l6l [2J]. The important, and so 
far unexplored, fact here is that (po varies significantly 
when the driving frequency is scanned across the value 
corresponding to the current reversal in the frequency 
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FIG. 1: Experimental results for ID rocking ratchet for cold 
atoms. The average atomic velocity, rescaled by the recoil ve- 
locity Vr {vr — 5.88 mm/s for *^Rb) is reported as a function 
of the frequency of the drive, for a fixed value of the rela- 
tive phase between harmonics of the ac drive. The different 
data sets correspond to different amplitudes of the driving 
force. The parameters of the optical lattice are: detuning 
from resonance A = — 9r and intensity per lattice beam: 
II = (43.5 ± 0.3) mW/cm^ The driving is a biharmonic 
force of the form of Eq. ([2]), with the relative phase kept 
fixed a.t (j) — TV /2 for the presented sets of measurements. The 
amplitude of the force is Fo = —mXgo where m is the atomic 
mass, A the laser field wavelength, and the values for go (in 
kHz'^) are, for the different set of data reported in the figures: 
(a) go = 3.2 • 10^ (b) go = 12.8 ■ 10^ (c) go = 19.2 • 10^ (d) 
go ~ 25.6 • 10^. The lines are a guide for the eye. 



domain. 

The dependence of the phase-lag (/jq and the amplitude 
A, as obtained by fitting data as those of Fig. [2] with the 
function v/vr = Asin{(j) — (j)o), on the driving frequency 
is reported in Fig. [31 The four different data sets corre- 
spond to the four different driving strength amphtudes 
for which the data in Fig. [l]were obtained. The phase- 
lag (jjQ shows a very large variation around the frequency 
at which the current reversal is observed in Fig. [1] with 
the phase lag varying from — tt to around the current 
reversal, and taking the value —tt/2 at the reversal fre- 
quency [27| . Since in the experiments reported in Fig. [1] 
the relative phase is fixed to (j> = 7r/2, this value of the 
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FIG. 2: Experimental results for ID rocking ratchet for 
cold atoms. The driving is a biharmonic force of the form 
of Eq. ([2j. The average atomic velocity, rescaled by the re- 
coil velocity Vr, is plotted as a function of the relative phase 
(j) between harmonics of the ac drive, for different values of 
the driving force frequency. The lines are the best fit of the 
data with the function v/vr ~ Asm{(ji — <^o). The lattice 
parameters are the same as in Fig. [T] The frequency of the 
drive for the different data sets is: (a) ijj/(2n) — 40 kHz, (b) 
uj/{2-k) = 50 kHz, (c) Lj/{2n) = 60 kHz, (d) uj/{2-k) = 70 
kHz, (e) t^/(27r) = 80 kHz, (f) uj/(2tt) = 100 kHz, (g) 
u)/{2-k) = 120 kHz. The driving force amplitude, for all data 
sets, is Fo — —mXgo with go ~ 19.2 ■ 10'^ kHz^. 



phase-lag, 0o — ~7r/2, guarantees that the current is 
suppressed at the reversal frequency, where the directed 
current inverts its sign. On the other hand, the am- 
plitude of the sin-like curve, see right panels in Fig. [31 
stays finite around the reversal frequency. Thus, we can 
conclude that the current reversal in the frequency do- 
main is determined by the large variation, around the 
reversal frequency, of the dissipation-induced symmetry- 
breaking phase lag (pQ. This generalizes the link between 
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FIG. 3: Left (a-d): dissipation-induced phase lag (po, as ob- 
tained by fitting data as those in Fig. [2l with the function 
v/vr = Asm{(j) — ^o), as a function of the frequency of the 
drive, for different values of the driving strength. Right (e-h): 
amplitude A, obtained from the same fits. Different rows of 
plots of the figure correspond to different driving strength, 
with the amplitudes of the drive the same as in Fig. [1] e.g. 
the data in (a) and (e) correspond to a driving amplitude 
go ~ 3.2 • 10'^ klfz'^ and so on. The lattice parameters are 
the same as in Fig. [T] The lines are a guide ot the eye. The 
arrows indicate the frequency at which the current reversal 
shown in Fig. [1] occurs. 



symmetry-breaking and current reversals established in 
Ref. ^6(1 in the case of current reversals in the amplitude 
domain. 



IV. NUMERICAL ANALYSIS 

The presented experimental results show that the cur- 
rent reversals in the frequency domain are determined by 
a sharp variation of the dissipative phase lag (f>o when the 
driving frequency is varied around the value correspond- 
ing to the current reversal. 

The conditions of our experiment correspond to the 
weakly damped regime, hence the appearance of a 
nonzero dissipation-induced phase lag (pQ. The valid- 
ity of the established link between current reversals and 
dissipation-induced symmetry breaking can be tested by 
considering the Hamiltonian and the overdamped limit. 
In both cases, the phase lag (po is fixed by the symmetries 
of the system valid in the respective limits: (j)Q = for 
the Hamiltonian limit and 0o — — 7r/2 [28j for the over- 
damped discussed in Sec. II. B. Thus (f>Q cannot 
vary with the driving frequency, and the current rever- 
sals observed in the case of moderate dissipation should 
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disappear. 

Our experimental set-up is not suitable to explore the 
two extreme case of no-dissipation, and very large dissi- 
pation. Thus, in order to show that these current rever- 
sals vanish when the phase lag (f)Q is fixed by the sym- 
metries in the Hamiltonian and overdamped regimes, we 
resort to numerical simulations of Eq. (1). 
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FIG. 4: (Color online) Average atomic velocity as a func- 
tion of the relative phase (p between harmonics of the ac 
drive for several values of the driving frequency uj, with 
ojy = k{2Uo/'m)^^^ and A\ = A2 = 1. Top panel: simulation 
data in the Hamiltonian regime (a = D = 0) for Fo — 0.2 Uok. 
Middle panel: simulation data in the weakly damped regime. 
The driving amplitude, the friction and the noise strength 
values are fixed to Fq = 0.2Uok, a — O.lSofo, and D = 
1.944_Do, respectively, where ao ~ mkvo, Do = UoVo/k, and 
vo = {Uo/my/^/10. The values about p = tt/2 (or p = 37r/2) 
show a current reversal as the driving frequency is increased. 
Bottom panel: simulation data in the overdamped regime for 
Fo = Uok, a = lOOofo, and D = 1.944 ■ 10^ Dq. Lines are a 
guide to the eye. 

Let us examine first the weakly damped regime with 
this model. The middle panel of Fig. 4 shows the aver- 
age particles' velocity as a function of the driving phase 
for different values of the driving frequency. A current 
reversal can be clearly observed at about p — n/2 (or 
<f) — 37r/2) as the driving frequency is increased, being 
a consequence of a variation of the phase lag 0o with 
the driving frequency. This is the regime explored in our 
experiment. 

On the other hand, the phase lag 0o is fixed to in the 
Hamiltonian regime due to the time-reversal symmetry. 
Accordingly, the top panel of Fig. 4 shows that no cur- 
rent reversals are observed when the driving frequency is 
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varied. In the limit of strong dissipation, the appearance 
of the supersymmetry fixes the phase lag to 0o = ~'''/2, 
as shown in the bottom panel of Fig. HI This fact pre- 
vents the occurrence of the current reversals observed in 
the presence of moderate dissipation. 

The numerical simulations thus confirm the validity of 
the established link between current reversals in the fre- 
quency domain and dissipation-induced symmetry break- 
ing. 

V. CONCLUSIONS 

In conclusion, we examined the mechanisms which de- 
termine current reversals in rocking ratchets as observed 
by varying the frequency of the drive. We considered the 
specific case of a bi-harmonically driven spatially sym- 
metric ratchet. We found that a class of these current 
reversals in the frequency domain are precisely deter- 
mined by dissipation-induced symmetry breaking. Cor- 
respondingly, these reversals are observed only for mod- 



erate dissipation, and disappear in the Hamiltonian and 
overdamped limit. 

Our experimental and theoretical work thus extends 
and generalizes the previously identified relationship fl^ 
between dynamical and symmetry-breaking mechanisms 
in the generation of current reversals. 

In the future, it would be interesting to re-examine the 
issue of the relationship between microscopic dynamics 
and current reversals in the system studied by Mateos 
[ist - Whether the current reversals discussed there can 
be associated to some distinctive features of the micro- 
scopic dynamics is still an open question. 
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